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Abstract 

The low-energy theorem for the forward Compton scattering is generalized to the 
case of an arbitrary target spin. The generalization is used to calculate the correspond- 
ing contribution to the deuterium hyperfine structure. The nuclear-structure correc- 
tions are quite essential in this case due to the deuteron large size. Corrections of this 
type calculated here remove the discrepancy between the theoretical and experimental 
values of the deuterium hyperfine splitting. Explicit analytical result is obtained also 
for the deuteron polarizability contribution to the Lamb shift. 
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1 Introduction 

The hyperfine (hf ) splitting in deuterium ground state has been measured with high accuracy. 
The most precise experimental result for it was obtained with an atomic deuterium maser 
and constitutes |T[j 

v exp = 327 384.352 522 2(17) kHz. (1) 
Meanwhile the theoretical calculation including higher order pure QED corrections gives 

vqed = 327 339.27(7) kHz. (2) 

The last number was obtained by using the theoretical result for the hydrogen hf splitting 
fromRef. @ 

1420 451.95(14) kHz 

which does not include proton structure and recoil radiative correction, and combining it 
with the theoretical ratio of the hf constants in hydrogen and deuterium from Ref. 

4.339 3876(8) 

based on the ratio of the nuclear magnetic moments and including the reduced mass effect 
in \ip(0)\ 2 . 

In the present article the discrepancy 

Vexp ~ vqed = 45 kHz (3) 

is removed by taking into account the effects due to the finite size of deuteron. Such effects 
are obviously much larger in deuterium than in hydrogen. One nuclear-structure contribution 
to the deuterium hf splitting was pointed out long ago ||] by some intuitive arguments, and 
then discussed in more detail in Refs. 0, |, 0. Here we treat the deuteron finite-size effects 
in a systematic way. Not only the old result is reproduced, but new corrections are obtained, 
among them that generated by the deuteron electric and magnetic form-factors. 

To calculate some contributions to the deuterium hf structure we generalize the low-energy 
theorem for the Compton scattering to an arbitrary target spin. 

One more subject considered in this paper is the contribution of the deuteron polarizability 
to the deuterium Lamb shift. The fact that this contribution is close to the accuracy attained 
in experiment was pointed out in Refs. [| where the effect was calculated in the square-well 
approximation for the nuclear potential. This correction was calculated then with separable 
nuclear potentials in Ref. flQfl . Here we obtain in the zero- range approximation a closed 
analytical result for the deuteron polarizability contribution to the Lamb shift. 



2 Low-energy theorem for the forward Compton scat- 
tering on an arbitrary target 

According to the well-known low-energy theorem for the Compton scattering on a spin 1/2 



hadron [11, 12], this amplitude is described by the pole Feynman diagrams. We are interested 
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here not in the spin-independent Thomson amplitude which is of zeroth order in the photon 
frequency uj, but in the next, spin-dependent term of the u expansion. This contribution 
as well can be easily obtained directly in the nonrelativistic approximation |13|| . In this 
approximation the electromagnetic vertex can be immediately written for an arbitrary spin 
s: 

e ( Z -> fi 



2m p {A 



(2p + k) + -i[s x k) 



(4) 



Here Z is the hadron charge, and its ^-factor is related as follows to the magnetic moment \i 
measured in the nuclear magnetons e/2m p : 

g = n/s. 

In the forward scattering case, when the hadron is at rest (p = 0) and both photons have 
physical, transverse polarizations ((ke) = (ke 1 ) = 0), this vertex reduces to the pure spin 
interaction. The nonrelativistic pole scattering amplitude generated by this interaction is 



Mi = M lmn e' 



2m r , 



ug 2 i(s- [e' x e]) 



(5) 



However, this expression is incomplete. Indeed, being applied to a proton, it does not repro- 
duce the well-known result [O, [12| according to which the spin-dependent forward scattering 



amplitude is proportional to (g — 2) 2 . The explanation was pointed out in Ref. ||13| : the 
nonrelativistic pole amplitude should be supplemented by a contact term generated by the 
spin-orbit interaction, which restores the agreement with the classical result |11, [12] . 

This contact term can be easily derived for an arbitrary spin (as well as the nonrelativistic 
pole contribution (|5|)) in the following way. The motion equation for spin in an external 
electric field E can be written to lowest nonvanishing order in v/c as 



ds 
~dt 



2m r 
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Z 
A 



sx [E x v\ 



(6) 



Here A is the target mass as measured in the units of m p (i.e., it is the atomic number in the 
case of nuclei we are mainly interested in in the present work.) The interaction Hamiltonian 
generating equation (0) is obviously 
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2m„ \ u A 



- (s-[Exv\). 



(7) 



Expressions (|6]), (0) are in fact only slightly rewritten formulae from book [14|. After substi- 
tuting 

p — ZeA 

v = : 

2Am p 

into Hamiltonian (^), we arrive at the following contact interaction: 

/ \ 2 
I p \ 

V, -- 



2m r 
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It produces one more piece in the scattering amlitude: 

M 2 = M 2mn e' m e n = u (-4 |) (g - | ) i (s- [e? x e]) . (9) 



Taken together, expressions (|5|) and (g) generate the forward scattering amplitude: 

M = (—) u (g - 2^Vi(s- [# x e]). (10) 



2m p J v A/ 



This is the generalization of the low-energy theorem we are looking for. 

In the particular case of a proton (s = 1/2, Z = A = 1) this formula reduces to the result 
ofRefs. P, pi. 



3 Low-energy theorem and deuterium hyperfine struc- 
ture 

Being dependent on nuclear spin, the low-energy amplitude obtained contributes to the 
atomic hyperfine structure. However, to apply it to this problem, the amplitude should 
be modified. Indeed, both photons exchanged between the nucleus and electron are off-mass- 
shell. So, here u ^ \k\. Then, virtual photons have extra polarizations. We will use the 
gauge A = where the photon propagator is 

d- k-k 

D im {UJ,k) = — =r-, di m = Si m — Dqq = Do m = 0. (11) 

ui — k u 
Now, first of all, the magnetic moment contribution M lmn to the pole diagram changes 

to 

M lmn = (77— ) 9 2 i£mnkk k (k ■ s) -■ (12) 
y 2m p J to 

Second, the convection current, which is proportional to ±/c for a nucleus at rest, is operative 
now and induces the following nuclear-spin-dependent contribution to the forward scattering 
amplitude: 

/ e \ 2 Z 1 

Mimn ( ~Z ] 2 — — Q i (k m € nrs k r S s /c n 6 mrs /c r S s ) . (13) 

V 2m J A 00 
We are ready now to write down the electron-nucleus nuclear-spin-dependent scattering 
amplitude generated by the two-photon exchange with the deuteron intermediate state: 

T el = 4vrm J P - 2lk ( lmn + 2mn + 3mn) ' ( ^ 

Here = (m e , 0,0,0) is the electron momentum. The structure — k + m e ) r yj reduces 
to —iujtijioi where a is the electron spin. We will calculate this Feynman integral with 
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logarithmic accuracy. Two subtle points should be mentioned here. The singularity at 
10 — 0, originating from 1/u 2 in the projection operator (|TTD, should be understood in the 
sense of the principal value. Then, one comes across a term containing integral 

d\k\ 



o \k\ 2 

which diverges in a power- like way at \k\ — > 0. To regularize it one should introduce non- 
vanishing electron velocity v ; in this way one obtains the well-known Coulomb wave-function 
correction ira/v which should be neglected to our purpose. In our logarithmic approximation 
one should neglect as well possible constant terms originating from this integral. 

The final result of calculation can be conveniently presented in the following form. Let 
us write down the spin-dependent Born term in the electron-nucleus scattering amplitude: 

T =-- g{o-s). 15 

6m e m p 

It is in fact minus Fourier-transform of the lowest order contact hyperfine interaction. There- 
fore, the ratio A e i = T e i/T is nothing else but the relative magnitude of the discussed 



correction to the hf structure. The result for integral flT4|) can be written as 

3a m e , A 1 / 2 Z Z 2 \ 

A el = — —log — - [g 2 -Ag- - 12— . 16 

87T m p m e g \ A A 1 J 

At s = 1/2, A = Z = 1 it agrees with the corresponding results [|15|, [16| for muonium 
(where the effective cut-off A is provided by the muon mass) and hydrogen (where the integral 
is cut off at the typical hadronic scale m p ). 

In the case of deuterium (s = 1, g — — 0.857, A = 2, Z = 1), we are mainly interested 
in, the integration over the momentum transfer k is cut off at the inverse deuteron size 
k = 45.7 MeV. In this way we obtain the following result for the relative correction in 
deuterium: 

A el = — —log— U d -2 . (17) 

87r m p m e \ n d J 

At larger momentum transfers, k > k, the amplitude of the Compton scattering on a 
deuteron is just the coherent sum of those amplitudes on free proton and neutron. This 
contribution to the hf structure can be also easily obtained with the above formulae. Since 
both nucleons are in the triplet state, one can substitute s/2 both for s p and for s n . With 
the logarithmic integral cut off here at the usual hadronic scale m p = 770 MeV, we get in 
this way 

3a m e , m n 1 . 9 . . 

A in = — ^ log^— //J- 2fi p - 3 + £). (18) 
Att m p k fid 

Here fi p = 2.79 and /i n = — 1.91 are the proton and neutron magnetic moments. 

In the conclusion of this section let us mention strong numerical cancellation between A e ; 
and A,-„.. 
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4 Contribution of deuteron virtual excitations 



The low-energy Compton amplitude discussed above is just linear term of the full amplitude 
expansion in u (and in \k\ 2 /u for virtual photons). One may expect, however, that for the 
deuteron with its small binding energy this approximation is insufficient even in the atomic 
problem considered here. And indeed, we will see below that the deuteron virtual excitations 
are far from being inessential to our problem, they strongly dominate the effect discussed. 
Since the contribution of large momentum transfers k > k has been calculated already (see 
formula (|TBP, we confine now to the region k < k. All calculations below are performed in 
the zero-range approximation (zra) for the deuteron which allows to obtain all the results in 
a closed analytical form. 

Let us start with the transitions induced by the spin interaction only. The corresponding 
scattering amplitude is 



M4 mn — — I ] ^ 



< 0\[k x S] m \n >< n\[k x S^] n \0 > 



2m p J V lu - E n - I 



< 0|[fc x S%\n >< n\[k x S\ m \0 > 
to + E n +1 



(19) 



Here / = k 2 jm p is the deuteron binding energy, E n = p 2 /m p is the energy of the intermediate 
state | n > (all intermediate states belong to the continuous spectrum), and 

S = fip<j p e ~ + fi n <j n e ^ 

where B p < n ) is the proton (neutron) spin operator. 

When calculating this contribution, we will retain only terms logarithmic in the parameter 
e = I / k = K,/m p 1. The log is generated by the integration over k, and to obtain it it 
is sufficient to put the exponents in S equal to unity. Then the operator S can induce only 
Ml transitions. In the zra the deuteron ground state is pure 3 Si from which Ml transition 
is possible also to S'-states only. But due to the orthogonality of the radial wave functions of 
different triplet states, the intermediate states confine to 1 S . 

Since the total spin operator s = (1/2) (c? p + a n ) does not induce triplet-singlet transi- 
tions, the operator S reduces here to 

In our problem of hf structure we need the antisymmetric part of tensor ( |19D which is linear 
in the deuteron spin s. It looks now as follows: 

mL = - fe) { ^ - emnkhCk ■ ^ u I ^ 2< - % p + wih m 

where < 1 5'o, p \ 3 Si > is the overlap integral of the ground state zra wave function 

r~K~ e~ Kr , , 

* = VS— (21) 
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with the singlet one of the momentum p. 

This contribution to the electron-deuteron scattering amplitude 



T 1 - Anrni f — dimd 3 n 7»(^ fc + m e )7j j 
m " 47Fm 7 (27r )4 fc 4 k 2 _ 2lk M mn W 

is easily calculated with the logarithmic accuracy. Indeed, to this accuracy the energy de- 
nominator in formula ( p0| ) can be simplified to 

1 



OJ 2 — K 4 /m 2 



Then the integration over p reduces to the completeness relation. The resulting relative 
correction to the deuterium hf structure is 

Al = 3a me mp (/i p - fi n ) 2 
8n nip k //,/ 

Let us consider at last the inelastic contribution induced by the combined action of the 
convection and spin currents. Since the convection current is spin-independent, all the in- 
termediate states are triplet ones as well as the ground state. Therefore, here operator S 
simplifies to 

S - s{^ 2 + ii n e-^' 2 ). (24) 

According to the common selection rules, neither 3 Si can be excited by the convection current 
from the ground state. But in the zra all the states with I ^ are free ones. It means that 
here we can use as the intermediate states just plane waves, eigenstates of the momentum 
operator. Then the only matrix element entering the amplitude is 



< \e ± ^' 2 \p> = . (25) 



(p ± k/2) 2 + k 2 
In this way the amplitude itself simplifies to 
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M±=\ — \ 2Kujr^-{- — ^ — — + 



dp j fi p H r , 



\2m p J ~ J it 2 { k/2) 2 + k 2 } 2 [{p- k/2) 2 + K 2 }[(p + k/2) 2 + k- 



(26) 



uj 2 — {p 2 + K 2 ) 2 /m 2 

Integrals we come across when calculating the corresponding part of the electron-deuteron 
scattering amplitude 

T 2 - W [ dmdjn 7i(f-fc + ™ e )7j M 2 (27 ) 

m ~ a J (2tt) 4 k 4 k 2 - 2lk m "' 1 j 
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are rather tedious even when we confine to terms singular in the parameter e = n/m p <C 1 : 
1/e and log e. This relative correction to the hf structure is to this approximation 

A 2 in = a — ^ — -log— ^ — . 28 

2k \i d n nip k /i d 

One of the terms in this correction, 

2k fi d 

was obtained and discussed in Refs. j|, [5], ^ [7|. However, another term 

m e n p 



a 

2k n d ' 



is larger numerically. 



5 Corrections due to finite distribution of the deuteron 
charge and magnetic moment 



In the case of hydrogen this problem was considered many years ago [I7|. For deuterium those 
corrections should be obviously larger. In the zra the problem has here a closed solution. 

Let us start with the second-order amplitude of the electron-deuteron scattering as in- 
duced by the deuteron charge and magnetic moment. The nucleus will be treated in the 
static limit. However, its finite charge and magnetic moment distributions will be taken into 
account by introducing the corresponding form-factors, F ch and F m . This amplitude is 

V = - {A.af^L f^-iisx q] 



2m p J (2tt) 3 l Hi q 4 

7o(j + q + m e ) 7 - 7 (1 + q + m e ) 70 
(I + q) 2 — ml 



(29) 



Here again = (m e , 0,0,0), and g M = (0,g). This expression can be conveniently trans- 
formed to 

8m e a f dq „ , 9 . „ , 9 . _ . . 

v = / -f F ch(q 2 ) F m (q 2 ) T (30) 

7r J q 

where T is the momentum- independent magnetic Born amplitude ([15]). 

The effect we are interested in, vanishes of course if unity is substituted for both form- 
factors. Therefore, the corresponding relative correction to the Born amplitude To and to 
the hf splitting is in fact 

8m P a f dq r 
A f = — ' 



7r J q 



J 4 [Fck(q 2 ) F m (q 2 ) - l] ■ (31) 
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In the zra both deuteron form-factors have simple form: 



F ch (q 2 ) = F m (q 2 ) = < 0|e^ 2 |0 > = j arctg£. (32) 

Substituting it into formula (]3Tf) , we get the following explicit expression for the correction 
discussed: 

Tfl 

A f = -a-i(l + 21og2). (33) 

Two closely related features of the effect (in no way specific for deuterium only) are 
worth emphasizing here. This correction is of first (but not second) order in the ratio of the 
nuclear size to the Bohr radius m e a/n. Then, contrary to possible naive expectations, the 
contributions of the charge and magnetic form-factors are not additive. Both circumstances 
can be traced back to the fact that the typical momenta entering integral (|3~I| ) are of the 
nuclear, but not atomic, scale. 

6 Deuterium hf structure, discussion of results 

Our total result for the nuclear-structure corrections to the deuterium hf structure, comprising 
all the contributions, (0), (0), fl2|), fl2|), is 

A= a^ {^-^ - -(l + 21og2)| + ^ m,^-^)* 
2 k [ Hd 3 j 8tt m p k fid 

3a m e nip fi p - fi n 3a m e k 1 2 

log — — + log (fj, d - 2fi d - 3 34 

7i nip k fid o7r nip m e fid 

3a m e nip 1 2 2 

+ log ( fi - 2 fa - 3 + fi n ). 

An nip k Hd 

Numerically this correction to the hf splitting in deuterium constitutes 

Au = 43 kHz. (35) 

It should be compared with the lacking 45 kHz (see @). Taking into account the approx- 
imations made, first of all the crude nuclear model (zra), then the neglect of nonlogarithmic 
contributions, we believe that the agreement is quite satisfactory. In particular, including the 
correction due to the effective interaction radius r into the normalization of the deuteron 
ground state wave function (see details in the next section) would certainly enhance some 
contributions. 

Clearly, the nuclear effects discussed are responsible for the bulk of the difference between 
the pure QED calculations and the experimental value of the deuterium hf splitting. The 
calculation of this hf correction, including accurate treatment of nuclear effects, would serve 
as one more sensitive check of detailed models of deuteron structure. 
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7 Nuclear polarizability and Lamb shift in deuterium 



The nuclear polarizability contribution to the Lamb shift was considered recently in Refs. 
H, [| Here we present an analytical calculation of the effect with a closed result. The 
zra approximation used by us is applicable when the region of the wave function localization 
is much larger than the interaction range. Essentially the same condition is necessary to use, 
instead of the true interaction, the crude approximation of the square-well potential, as it is 
done in Refs. g, ||. 

The effect we are interested in now, is due to the photon-deuteron scattering amplitude 
induced by the convection current only. We will see that in this problem, that of the nu- 
clear polarizability contribution to the Lamb shift, the characteristic values of the photon 
4-momenta are as follows: 

\k\ </= — \p\. (36) 

m p 

Therefore, now we can omit in the Compton amplitude all dependence on k. As well as in the 
amplitude M^ n , all the intermediate states here have I ^ and can be described therefore by 
plane waves. We will use again matrix element (p5|), but this time at k — 0. At last, in the 
present problem we are interested in the scalar part of scattering amplitude which reduces to 

P 2 f 1 1 1 




\2m p J 3 J 7r 2 (p 2 + k 2 ) 2 \lo — (p 2 + K 2 )/m p u + (p 2 + K 2 )/m p 

We subtract from the expression in braces the term 



-2 mp 



p 2 + K 



After integrating over p this term being added to the Thomson scattering (seagull) amplitude 
for a proton, — e 2 /m p , reproduces the correct one for a deuteron, — e 2 /2m p . With the identity 

1 1 2 2u 2 



CO — It UJ + u u (u 2 — u 2 ) u 
we get the following expression for the photon-deuteron scattering amplitude in question: 



(37) 



7T 2 (p 2 + K 2 ) 3 {uJ 2 — (p 2 + K, 2 ) 2 /m 2 } 

The contribution of this tensor to the electron-deuteron scattering amplitude 

T 3 - d.irm ( d * k dimd i n JiV - fc + m e )7j 3 . . 

m - m J (2vr) 4 k A k 2 - 2lk mn [ } 
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can be easily transformed to 



T] 



32-K 2 a 2 K f dp p 



;2 



3m p J vr 2 (p 2 + n 2 f 



d 4 k [ I 2lu 3 1 1 

l.A 1 1,9 n/L\ ( ~ 9 7~3. \ ..9.\9. I 9 ' V**) 



(2tt) 4 \u(k 2 - 2lk) k 4 (k 2 - 2lk)j lu 2 - (p 2 + K 2 ) 2 /m 2 p 

The first term in the braces here contains no photon propagation, neither l/k A , nor 1/k 2 . 
In other words, it corresponds to the instantaneous Coulomb interaction. The second term 
corresponds to the exchange by three-dimensionally transverse quanta, i.e., to the magnetic 
interaction of convection currents. 

Perhaps, the most convenient succession of integrating expression (R^) is as follows: the 
Wick rotation; transforming the integral over the Euclidean u to the interval (0, oo); the 
substitution k — > ku; integration over u; integration over k (at the last two procedures it 
gets clear that the effective values of u>, \k\ belong to interval (|36D indeed); at last, integration 
over p. The following identity is useful here: 

'^(l-xrv-Mogx = ^ a) ^ ( ,? [m - ^(a+b)}-, m = llogr(6). 

o L (a + b) do 

The effective electron-nucleus interaction operator (equal to — Tf n ) can be finally pre- 
sented in the coordinate representation as 

V le = -«m e a d (0)5flog— + ^]6{r). (40) 
V m e 20/ 

Here a^(0) is the static value of the deuteron electric polarizability defined as usual by the 
relation 

2 r dp p 2 + k 2 < 0\r \n >< n\f |0 > 
a d (uj) = Ana - J — — — ^ - _ -. (41) 

The matrix elements entering expression ([|l]) are dominated by large distances. In this 
asymptotic region the naive zra expression for the deuteron ground state wave function 
should be supplied by the correction factor (1 — ro^)^ 2 taking into account finite effective 



interaction radius r (see Refs. |L4], [T^j). In this way we get the following result for the static 
electric polarizability: 



This numerical value is close to the experimental one |T3: 0.70(5) fm 3 (as well as to the values 
0.613, 0.623, 0.625 fm 3 obtained in Ref. [|K]] with different separable nuclear potentials and 
to 0.635 fm 3 found in Ref. |] with a square- well potential). 

The overall result ( pEUD consists of two contributions of different physical origin. The 
dominating one is generated by the instantaneous Coulomb interaction. Its contribution to 
the overall numerical factor 

- 5 l l0g ^ + 20 
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in formula (HO) is 

/ 8/ 
-4 log— + - 

V m e 12, 

The magnetic interaction contributes to the overall factor 

, 8/ 17 
log — 

m e 12 

The level shift of the deuterium ground state produced by operator ([10) constitutes 

— 22.3 kHz. The Coulomb and magnetic contributions to it are, respectively, — 19.7 and 

— 2.6 kHz. The results are close to those of Refs. [§, |, 0- 

No wonder that the Coulomb contribution is negative: this is a true second-order (in the 
electron-nucleus static interaction) correction to the ground state of a system consisting of 
an electron at rest and a nucleus which is in the ground state itself. The sign of the mag- 
netic contribution cannot be fixed in this way: in the language of the common noncovariant 
perturbation theory this is a fourth-order correction, second-order in the photon-electron 
interaction and second-order in the photon-nucleus one. 

One more contribution to the Lamb shift in deuterium is caused by the deuteron magnetic 
polarizability, considered earlier also in Ref. |I0|| . This is in fact the contribution of the 
scalar part of amplitude (|19|). The calculation simplifies due to the following circumstances. 
First, the numerators di m of the photon propagators reduce here obviously to 5i m . Then, 
the integration over k is spherically-symmetric one. So, to our purpose the scalar part of 



amplitude ( |19D may be simplified to 

(£*p - fi n ) 2 k(k + ^i) 2 



M = - Ana 

mn 9ml 



p 



dp 



7T 



2 (p 2 + k, 2 ) (p 2 + k. 2 ) [uj 2 — {p 2 + K 2 ) 2 /m 2 



We have used here the explicit form of the x Sq coordinate wave function in the zra: 

V 2 ixr 

where 

ctg5 = — , ki = 7.9MeV. 
p 

Its overlap with the ground state zra wave function ( pi] ) constitutes 



<'S f Sl >= ^iH±M . (45) 

(p 2 + K 2 )Jp 2 + k\ 



Further calculations are close to those related to the electric polarizability; only the 
last integration, that over p, is done numerically for the nonlogarithmic contribution. The 
resulting effective electron-nucleus interaction operator can be written as 

V lm = am e (3 d (0) flog— - 1.24) 5(r). (46) 
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Here /3<f(0) is the static value of the deuteron magnetic polarizability equal to 



&(0) 



a {jJL p - fi n ) 2 1 + Kx/3/C 



i m p k z 



(47) 



1 + K\/k 

This contribution to the Lamb shift of the deuterium ground state constitutes 0.31 kHz which 



is very close to the result of Ref. 10 
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